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REACTIONS 
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Abstract. Many phenomena in biology involve both reactions and chemotaxis. These pro- 
cesses can clearly influence each other, and chemotaxis can play an important role in sustaining 
and speeding up the reaction. However, to the best of our knowledge, the question of reac- 
tion enhancement by chemotaxis has not yet received extensive treatment either analytically 
or numerically. We consider a model with a single density function involving diffusion, advec- 
tion, chemotaxis, and absorbing reaction. The model is motivated, in particular, by studies 
of coral broadcast spawning, where experimental observations of the efficiency of fertilization 
rates significantly exceed the data obtained from numerical models that do not take chemotaxis 
(attraction of sperm gametes by a chemical secreted by egg gametes) into account. We prove 
that in the framework of our model, chemotaxis plays a crucial role. There is a rigid limit to 
how much the fertilization efficiency can be enhanced if there is no chemotaxis but only advec- 
tion and diffusion. On the other hand, when chemotaxis is present, the fertilization rate can 
be arbitrarily close to being complete provided that the chemotactic attraction is sufficiently 
strong. Moreover, an interesting feature of the estimates on fertilization rate and timescales in 
the chemotactic case is that they do not depend on the amplitude of the reaction term. 



1. Introduction 

Our goal in this paper is to study the effect chemotactic attraction may have on reproduction 
processes in biology. A particular motivation for this study comes from the phenomenon of 
broadcast spawning. Broadcast spawning is a fertilization strategy used by various benthic 
invertebrates (sea urchins, anemones, corals) whereby males and females release sperm and egg 
gametes into the surrounding flow. The gametes are positively buoyant, and rise to the surface 
of the ocean. The sperm and egg are initially separated by the ambient water, and effective 
mixing is necessary for successful fertilization. The fertilized gametes form larva, which is 
negatively buoyant and tries to attach to the bottom of the ocean floor to start a new colony. 
For the coral spawning problem, field measurements of the fertilization rates are rarely below 
5%, and are often as high as 90% [HI [151 [26l [31] • On the other hand, numerical simulations based 
on the turbulent eddy diffusivity [1] predict fertilization rates of less than 1% due to the strong 
dilution of gametes. The turbulent eddy diffusivity approach involves two scalars that react 
and diffuse with the effective diffusivity taking the presence of the flow into account. It is well 
known, however, that the geometric structure of the fluid flow lost in the turbulent diffusivity 
approach can be important for improving the reaction rate (in the physical and engineering 
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literature see |2H |28l [33]; in the mathematical literature see [HI |3l [131 [S [E] for further 
references). Recent work of Crimaldi, Hartford, Cadwell and Weiss [HI [Z] employed a more 
sophisticated model, taking into account the instantaneous details of the advective transport 
not captured by the eddy diffusivity approach. These papers showed that vortex stirring can 
generally enhance the reaction rate, perhaps accounting for some of the discrepancy between 
the numerical simulations and experiment. 

However, there is also experimental evidence that chemotaxis plays a role in coral fertilization: 
eggs release a chemical that attracts sperm [H [21 [181 [19] . Mathematically, chemotaxis has been 
extensively studied in the context of modeling mold and bacterial colonies. Since the original 
work of Patlak [2S] and Keller-Segel [HI [T2] where the first PDE model of chemotaxis was 
introduced, there has been an enormous amount of effort devoted to the possible blow up and 
regularity of solutions, as well as the asymptotic behavior and other properties (see [27] for 
further references). However, we are not aware of any rigorous or even computational work on 
the effects of chemotaxis for improved efficiency of biological reactions. 

In this paper, we take the first step towards systematical study of this phenomenon, by 
analyzing rigorously a single partial differential equation modeling the fertilization process: 



Here, in the simplest approximation, we consider just one density, p(x, t) > 0, corresponding 
to the assumption that the densities of sperm and egg gametes are identical. The vector field u 
in (11. ip models the ambient ocean fiow, is divergence free, regular and prescribed, independent 
of p. The second term on the right is the standard chemotactic term, in the same form as it 
appears in the (simplified) Keller-Segel equation (see [27j). This term describes the tendency 
of p{x, t) to move along the gradient of the chemical whose distribution is equal to —/Sr^ p. 
This is an approximation to the full Keller-Segel system based on the assumption of chemical 
diffusion being much faster than diffusion of gamete densities. The term (— ep^) models the 
reaction (fertilization). The parameter e regulates the strength of the fertilization process. The 
value of e is small due to the fact that an egg gets fertilized only if a sperm attaches to a certain 
limited area on its surface (estimated to be about 1% of the total egg surface in, for example, 
sea urchins eggs [30]). We do not account for the product of the reaction - fertilized eggs - 
which drop out of the process. We are interested in the behavior of 



which is the total fraction of the unfertilized eggs by time t. It is easy to see that mQ{t) is 
monotone decreasing. High efficiency fertilization corresponds to mo(t) becoming small with 
time, as almost all egg gametes are fertilized. We prove the following results. 

Theorem 1.1. Let p{x,t) solve (11. ip with a divergence free u{x,t) G C°°(M'^ x [0, oo)) and 
initial data po > G iS(]R^) (the Schwartz class). Assurae that qd > d + 2, and the chemotaxis 
is absent: % = 0. Then there exists a constant po depending only on e, q, d and po{x) but not 
on u{x, t) such that mo(t) > po for allt>0. 

Moreover, po — )■ mo(0) as e — )■ while po,u and q are fixed. 



hp + u-Vp = Ap + xV(pV(A)-V) - ep^ pi^, 0) = po{x), x E M^. 



(1.1) 




BIOMIXING BY CHEMOTAXIS AND ENHANCEMENT OF BIOLOGICAL REACTIONS 3 

Remarks. 1. Observe that the constant fiQ does not depend on u. No matter how strong the 
flow is or how it varies in time and space, it cannot enhance the reaction rate beyond a certain 
definitive hmit. 

2. The condition qd > d + 2 does not include the most natural case of d = q = 2. Dimension two 
corresponds to the surface of the ocean, and q = 2 corresponds to the product of egg and sperm 
densities. Our preliminary calculations show, however, that the mathematics oi d = q = 2 case 
is different and more subtle. Then the norm of p for sufficiently rapidly decaying initial 
data does go to zero but only very slowly in time. The difference between chemotactic and 
chemotactic-free equation (II. ip in this case is likely to manifest itself in the time scales of the 
fertilization process: in the presence of chemotaxis the norm initially decays much faster. 
We will address this issue in a separate publication, to keep the present paper as transparent 
as possible. 

3. The case of small e is interesting due to experimental relevance. It will be clear from the 
proof that decrease of the norm, mo(0) — limf_j,oo ^o(^) is of the order e when chemotaxis is 
absent. 

4. The condition that po & S can of course be weakened. What we need is the initial data 
that is decaying sufficiently quickly and is minimally regular. Similarly, the condition that u is 
smooth can be weakened to, say, without much difficulty. We decided not to pursue most 
optimal regularity conditions on po and u in this paper to simplify presentation. 

5. By u G C°°(]R X [0, oo)) we mean that bounds on every derivative of u are uniform over all 
(x, t) eR'^ X [0, T], for every T > 0. 

On the other hand, in the presence of chemotaxis, we have 

Theorem 1.2. Let p{x,t) solve (11. ip with a divergence free u{x,t) G C°°{R'^x [0,oo)) and fixed 
initial data po > E S. Assume that d = 2, and q is a positive integer greater than 2. Then we 
have that mo{t) — )■ c(x,Po5w) > as t —)■ oo, hut c(x,Po5^) — )■ as x — ^ C)0, with q, po and u 
fixed. 

Remarks. 1. In general, chemotaxis can lead to finite time blow up of the solution (see 
[27] for references). However, it is known that the presence of the reaction can lead to global 
regularity (see [32] for a slightly different model where references to some earlier works can 
also be found). We will sketch the global regularity proof for solutions of the equation (II. ip in 
Appendix I. 

2. We prove more (see Theorem 14. 2p . Here we stated the result in the simplest form to avoid 
technicalities. 

3. An interesting feature of Theorem 11.21 is that c{x,Po,u) — ?■ uniformly in reaction strength 
e. The timescale to achieve most of decay also does not depend on e (Theorem 14. 2p . The 
mechanism of this effect will be clear from the proof. Roughly speaking, what happens is the 
solution tries to blow up to a 5 function profile. The reaction term, however, prevents blow 
up from happening. However, the solution gets larger if e is small. Chemotaxis and reaction 
balance in a way to produce same order reaction effect independently of e. 

4. The case c? > 2 is mathematically different and it is not clear that the norm may become 
arbitrarily small in this case even with strong chemotaxis aid. There appears to be a genuine 
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mathematical reason why coral gametes rise to the surface instead of trying to find each other 
in the three-dimensional ocean! 

Hence our model implies that the chemotactic term, as opposed to the flow and diffusion 
alone, can account for highly efficient fertilization rates that are observed in nature. Moreover, 
Theorems 11.11 and 11.21 suggest that the presence of chemotaxis may be a necessary and crucial 
aspect of the fertilization process. Of course, a more realistic model of the process is a system 
of equations involving two different densities. We will show that even for the system case, the 
flow can only have a limited effect on fertilization efficiency, similarly to our simple model. It 
is possible that in the system case the flow and chemotaxis can play supplementary role, with 
flow acting on larger and chemotaxis on smaller length scales. Note that one can expect that 
in the system setting the chemotaxis effect is weaker, since it only appears in the equation 
for the sperm density. The influence of chemotaxis in the system setting, and investigation of 
quadratic reaction term are left for a later study. 



2. The reaction-advection-diffusion case 
In this section, we prove Theorem 11.11 Consider equation fll.ip with x = : 

9tp + M- Vp = Ap-ep^ p(x,0) =po(x). (2.1) 
As the first step, observe that by comparison principle, p{x,t) < b{x,t), where 

dtb + u-Vb = Ab, b{x,0) = po{x). (2.2) 
Also, note that since p{x,t) > 0, 



dt\\pi-,t)\\L^ = dt / p{x,t)dx = -e / p''{x,t) dx > -e I 

jRd J^o 



¥{x, t) dx. 



Therefore, the behavior of the L'^ norm of b can be used for estimating decay of the norm 
of p. We have the following lemma, similar in spirit (and proof) to Lemma 3.1 of [9]. 

Lemma 2.1. There exists C = C{d) that, in particular, does not depend on the flow u, such 
that 

||6(-,t)|U2 < min(||6o||L^Cr'^/1l&o||Li), \\K-,t)U- < min(||6o(x)|Uo., Cr^/2||6o||^,). (2.3) 
Proof. By Nash inequality [2T], we have 



\\bC-^<c{d)\\b\\%'\\vbU.. 

Multiplying f l2.2p by b, integrating, and using incompressibility of u, we get 

2 \\b\\h WboWh 

We used the conservation of the L^-norm of b in the last step. Set z{t) = \\b{-,t)\\j^2- Then 

z\t)<-cz{ty+^\bofJ. 
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Solving this differential inequality, we get 



z{t)< 



2Ct 1 

+ 



-d/2 



d\\ho\\li ||6o|ii2 
implying 

||6(-,t)||i.<min(||6o||i.,C7(rfr'^/^||6o|liO- 
This gives the first inequality in fl2.3p . 

The second inequality in (12 .Sp follows from a simple duality argument using incompressibility 
of u. Indeed, consider 6{x, s), a solution of 

dsO + u{x, t-s)-Ve = AO, e{x, 0) = Ooix) G S. 

A direct calculation shows that 

d f 

— / b{x, s)6{x,t — s) dx = 0. 
as Jj^d 

When s = t, we get 

b{x, t)6Q{x) dx 
For s = 0, this implies 

bQ{x)9{x, t) dx 



(2.4) 



< \\bix,mL4do\\L^ < cid)t-^/^bo\\LAML^- 



< c{d)t-''/4bohAML^ 

for every 60,^0 ^ S. Hence 

mx,t)\\L^<c{d)t-''/4eo\\L^ 

for every Oq E L^. To finish the proof of the Lemma, given t > 0, note that 

mx,t)\\L- < C(d)(t/2)-'^/^||6(x,t/2)|U. < C(rf)r^/2||6o||Li. 
Here in the second step we used (12 .40 and adjusted C{d). 



□ 



For a more precise estimate on the residual mass /iq, we need one more lemma. 

Lemma 2.2. Assume that p{x,t) solves ()2.ip with a smooth bounded incompressible u and 
Po G S. Then for every t > we have 

\\pix,t)\\LP / WPoWlp 



< 



for all 1 < p < 00. 



\\p{x,t)\\Li llPollLi 

Proof. For p = 1 the result is immediate. Consider some 1 < p < 00, and look at 



d ( j^dp^dx 



pi pdx 



-p-i 



/ p^'^{-u-Vp + Ap- tp'^)dx I pdx- I pPdx {-u ■ V p + Ap - ep'^) dx 
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Consider the term in the second hne above, which after integration by parts simphfies to 

(-{p-1) [ ff-^\Vp\'^dx-e [ p'^+P-Ux] [ pdx + e [ p^ dx [ p" dx. 
This does not exceed 

-e / p^+P-^dx I pdx + e I ff dx f p'^dx, 

which is less than or equal to zero by an application of Holder's inequality. 

The p = oc case follows by a limiting procedure since p{x, t) E S for all t. □ 

We are ready to prove Theorem 11.11 

Proof of Theorem The idea of the proof is very simple. We will show that if L^-norm of 
p at some time to is sufficiently small then for all times t > to the L^-norm of p{x,t) can not 
drop below ||p(to)IUi/2. This shows that p{x,t) can not tend to zero as t — )■ +00. 
Recall that for every t, 



dt / p{x,t)dx = —€ / p{x,ty dx > —e / b{x,ty dx, 

JR'i JR'i JR'i 

where b is given by (12. 2p . By Lemma [2. II and Holder's inequality, 

/ b{x,tydx < Cmin (\\po\\%,t-'-^\\po\\l) . 



Thus, for every r > 0, 



dt \ b{x,tydx<C{d)\\\p{-M)VL.r+\\p{-M)\\l. / {t-toY^dt 

to JR'i V Jto+T 

<C(rf,g)(||p(-,to)rL;o^||p(-,to)||L.r+||p(-,to)|ll.r^). (2.5) 

We used the assumption qd > d + 2 when evaluating integral in time. 

Assume, on the contrary, that the norm of p does go to zero for some u. Consider some 
time to > when ||p(-, to)||Li is sufficiently small (we'll have a precise bound later). Using 
Lemma [2.21 and (12. 5p . we see that further decrease of the norm from that level is bounded 

by 

||p(-,to)|Ui - ||p(-,t)|Ui < C{d,q)e (^^||p(-,to)|li.r+ ||p(-,to)||i.r^] , (2.6) 
for all t > to, T > 0. Choosing r to minimize the expression ( 12. 6p . we find that for every t > to, 

qd-d-2 

l|p(-,to)||Li-||p(-,t)IUi<C7(g,rf)e||p(-,to)|llifM^^l ' • (2-7) 
If ||p(-,t)||ii — )• as t — ^ +00, we may choose to so that 



qd-d-2 

C(g,ci)e||p(-,to)|ll7MT4r^) ' A' ^^.8) 
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Then we get that 

1 / 1 If IIpoIIli \ ^^^^^ 

\\p{-,t)\\L^ > 7:||p(-,to)||Li > Ho{q,d,po) = mill -||po||li 



2lirv 7"u/i|j^ — r'^yj \ ""7 r / i r\ wr'o : g i i i 

29-1 e-j-i C(g, d) 9-1 V||po||L°= 

for every t > to- This is a contradiction to the assumption that ||p(t)||Li — )■ as t — +oo. This 
argument can also be used to define po in the statement of the theorem. The last statement 
of the theorem is easy to prove by changing the condition < ^ in fl2.8p to < k where k can be 
taken as small as desired. □ 



3. The reaction-advection-diffusion case: a system 



In this section we show that the results of Section [2] largely extend to a more general model. 
Consider the following system 

dts = {u-V)s + KiAs ~ e{sey^^, s{x,0) = so{x) (3.1) 
a^e = (u ■ V)e + Ae - e(se)^/^ e(a;, 0) = eo(a;). (3.2) 

Here s{x,t) and e{x,t) are sperm and egg densities respectively. The following analog of The- 
orem 11.11 holds. 

Theorem 3.1. Let s{x,t), e{x,t) solve fl3.ip . fl3.2l) with divergence free u{x,t) E C°°(M'^ x [0, oo)) 
and initial data Sq, Cq G S . Assume that qd > d + 2, q > 2 and the chemotaxis is absent: x = 0. 
Then there exists a constant /ii depending only on e, q, d and eo(a;), So(a;) such that the 
norms of s{x,t) and e{x,t) remain greater than pi for all times. 



Remarks. 1. The condition q > 2 can be omitted if ||so||li 
2. Similarly to Theorem I l.lj one can show that limt_j.oo \\s{-, t) \ 
II Co II Li provided that q, u, Sq and Cq remain fixed. 



IIpoIIli- 



Li 



soIIli and limt_^oo ||e(-,t)||Li 



Proof. As before, we know that s{x,t) < s{x,t) and e{x,t) < e{x,t) where s,e solve (12.21) with 
initial data sq and eo, and the diffusion coefficients ni and K2, respectively. Lemma |2TT] can still 
be used to control s, e. Instead of Lemma 12.21 we will use a cruder bound. 

Observe that if ||eo||Li 7^ II-SoIIlI; then the norm that is larger initially remains larger than 
the other norm. Hence, assume without loss of generality that HcoUli < ||so||li and focus on 
the decay of ||e(-, t)||ii. Let us estimate the decay after some time to '■ 



dt 



to 



tQ+T 

to 

< t\\s( 



dt 



s{x,ty/^e{x,ty/^ dx < 

s{x,ty/^e{x,ty/^ dx 

^o)||fc'||e(-,to)IUi 
ei;to)\\L^+CT'-'- 





poo r 


+ 


I dt I 




Jto+T Jw 



19/2 



s{x,ty/^e{x,ty/^dx 

\s{;t)ri!M-Mi"dt 



n^-l 

2 

I Co 1 1 Loo 



to+T 

^ Ipollii 



(3.3) 
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Choosing r to minimize (13 .Sp leads to 

q{qd-d-2) q (q-2)(qd-d-2) , q-2 

\\e{-M)\W - U-Ml^ <C{q,d)e\\s4^'^'^-'' ||so||^r^1eo|L J"''"^' ||e(-, to) L/''^^ (3.4) 

Suppose that ||e(-,t)||Li does go to zero as t — t- oo. Choose to so that C||e(x, to)!!]^!'"^^ < \ 

1+ 1-^ 

(where C is the constant in front of ||e(-,to)||^i ^' in fl3.4p ). In this case, due to (13 ■4p . the L} 
of e{x,t) can never drop below half of its value at to- This is a contradiction. □ 

4. Reaction enhancement by chemotaxis 

In this section, we will show that chemotaxis, as opposed to a divergence free fluid flow, can, 
in principle, make reaction as efficient as needed. We consider the equation 

dtp = Ap-{u- V)p + xV(pV(A)-V) - ep^ p{x. 0) = po(x). (4.1) 

We will prove that the large time limit of the norm of p{x,t) goes to zero as chemotaxis 
coupling increases, independently of e. On the other hand, we will also prove lower bounds 
showing that the norm does not go to zero as t — )■ cxd for each fixed coupling. Before we 
state the main results of this section, there is an auxiliary issue we need to settle. In general, 
solutions to the chemotaxis equation may lose regularity in a finite time (see e.g. |27] for further 
references). As Theorem 14.11 below shows, this does not happen with the additional negative 
reaction term — ep*?, g > 2 in the right hand side: solutions with smooth initial data stay 
smooth. We will work with initial data which is concentrated in a finite region, in particular, 
with a finite second moment. As we will see, this property is also preserved by the evolution. 
Let us define 

ii/iim„= / (iv/i + i/(x)i)(i + |xr)dx. 

Let H"^ denote the standard Sobolev spaces in R''. Define a Banach space Kg n with the norm 
ll/lk.„ = ll/b» + ll/IUf„. Then we have 

Theorem 4.1. Assume that q > 2, n > and s > d/2 + 1 are integers and po G Ks,n- Suppose 
that u G C°°(]R'^ X [0, oo)) is divergence free. Then there exists a unique solution p{x,t) of the 
equation dH]) in C{Ks,n, [0, oo)) n C°°(R"' x (0, oo)). 

The proof of Theorem 14.11 uses fairly standard techniques; we sketch it in Appendix I. 
First, we prove the bound showing reaction enhancement by chemotaxis. Let us define 

m2 = Toaiuxg / \x — xo\'^poix) dx. 

Theorem 4.2. Let d = 2, and suppose that u G C°°(]R'^ x [0,oo)) is divergence free. Assume 
that q > 2, s > d/2 + 1 and n > 2 are integers and p{x, t) solves (14. ip with po > G Kg n- Then 
a. If u = 0, then lim(_^oo ||p(") i^) IIl^ < '^X"^- More precisely, for every r > 0, we have 

||p(,,)||,.<^(l + yr^. (4.2) 
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b. If u 7^ 0, then limt_j.oo ||p('5^)||li < C{u,m2)x '^^^- Moreover, for < r < x^^^ we have 

M;r)\\L.<Ciu,m,)ixr)-'/'. (4.3) 

Remark. Note, in particular, that if u = 0, the level ||p(-, t)!!/,! ~ x~^ will be reached in at 
most T ~ X, while the level ~ iii most r ~ 1. If m 7^ 0, the upper bound on the time 

scale to reach the norm level ~ is also r ~ 1. 

Proof. Since po ^ -^s.n, there exists such that f^^ \x — Xo\^po{x) dx = m2. Set = for 
simplicity. Consider 

dt I \x\^pdx= I \x\^{u-V)pdx + I \x\^^pdx + x / \x\^V {pV p) dx — e I \x\^p'^dx. 

(4.4) 

Observe that due to V ■ m = 0, 



/ |xp(M ■ V)pdx = -2 / 



■ u)pdx, 



and in dimension two 

/ |xpApda; = 4 / pdx. 

JR2 

For the chemotaxis term, we have 



/ \x\'^V{pVA-^p)dx = -2 [ ^ "^^ ^\ {x,t)p{y,t)dxdy = - ( [ pdx^ . 

iR2 JR2xM2 \X — y\ \Jj^2 J 

In the last step, we used symmetrization in x,y. Due to Theorem 14. all integrations by parts 
are justified for all t > 0. Therefore, we can recast f l4.4p as 

dt / Ixl"^ pdx = —2 / {x-u)pdx + A / pdx — x[ pdx\ — e / \x\'^p'^dx. (4.5) 

iR2 JR2 Jr2 \Jr'2 J iR2 

First let us set m = in fj45|) . Suppose that \\p{-,t)\\Li > Y for all t G [0,r], and F > 4/x. 
It follows from (14.51) that we need tY{x^ — 4) < m2 to avoid contradiction that m2 vanishes. 
This quadratic inequality translates into (14. 2p . 

Now, assume that u is an arbitrary smooth divergence free vector field. In this case, we 
further estimate 

with /3 > to be chosen. Then, it follows from (14. 5 p that 

dt / \x\'^pdx < 2x^^ / \x\'^ p dx + i A + 2x^ \\u\\\oc — X / pdx) pdx, 

and thus 

Stfe-'^""*/ \x\^pdx) Ke-^""'^' (4 + 2x^\\u\\l^-x [ pdx]! pdx. (4.6) 

\ iR2 / \ iM2 / Jr2 
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Assume now that for all t G [0,r], we have ||p(-,t)||ii >Y>0, and that 

Y>-{2 + x^\\u\\l^). 
X 

Then, the integral in time of the right hand side in (14. 6p over [0, r] can be estimated from above 
by 

l\-'^-''Y{4 + 2x^\\u\\l^-xY) dt={l-e-'^-'-)Yx^{2 + x^\\u\\U-xY/2). (4.7) 

Setting r = x^, we see that to avoid a contradiction, we need 

(1 - e-2)/r {xY-4- 2x''\\u\\l^) < 2m2. (4.8) 

An elementary computation shows that the optimal choice that makes Y the smallest is /3 = 1/3. 
Solving this quadratic inequality, we find that ||p(-,r = X"'^''^)IUi cannot exceed c{u,m2)x~'^^^ ■ 
More generally, for < r < x^''^, we get from 04.71) the bound 

||p(-,r)|Ui <C(w,m2)(xr)-i 

□ 

Observe that the reaction term did not play any quantitative role in the estimates. In 
particular, all estimates of the norm decrease and timescales are independent of e. The only 
role the reaction term plays is making sure we have a smooth decaying solution, so that all 
integrations by parts are justified. On a qualitative level, what happens is chemotaxis without 
reaction would lead to a 5 function profile blow up. With reaction present, the growth in the 
L°° norm of the solution is determined by the balance between chemotaxis and reaction term. 
Weaker fertilization coupling parameter e leads to stronger aggregation due to chemotaxis and 
thus effectively the same fertilization rates. In particular, the amount of the density that does 
react satisfies the lower bound independent of e. The same holds true for the time scale on 
which this reaction takes place - it is bounded from above by the pure chemotaxis blow up time 
independently of e. 

Next we prove a result in the opposite direction, showing that at least some estimates of 
Theorem 14.21 scale sharply in x- 

Theorem 4.3. Let d = 2, and suppose that u G C°°(]R'^ x [0,oo)) is divergence free. Assume 
that q > 2, s > d/2 + 1 and n > 2 are integers and p{x, t) solves (14. ip with po > G Kg^n- Then 
limt_^oo ||p(") i^) > 0. Moreover, for some initial data po, \\p{-,t)\\Li remains above c{q,po)x~^ 
for all times. 

Proof. Recall that dt J^2 p{x,t) dx = — J^2 p{x,tY dx. Let us derive estimates on HpUl"?- Multi- 
plying (14. ip by p'^~^ and integrating, we obtain 



Observe that 



-dt [ pUx= [ p^-^Apdx + x [ p''~^V ■ip\/A~^p)dx-e f p^^'Ux. (4.9) 
"'^V ■{pVA-^p)dx = -iq-l) [ p'''^Vp-VA-^pdx='^^ f p'^+Ux. 
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The last equality is obtained by integration by parts. Thus, we can rewrite (14. 9p as 

dt [ p^dx = / \^pi/^\^dx + x{q-l) [ p^^^dx-qef p^^^Ux. (4.10) 

Let us introduce v = p'^/^, and recall a Gagliardo-Nirenberg type inequality 

\\v\\l.^. < Cid,a)\\Vv\\^\\v\\^, (4.11) 

which is valid for all a > 0, > 1. In our case, we set a = 2/g, and inequality (14.111) translates 
into 



/ p''+^dx<Ciq) [ \Vp''^^\^dx [ pdx. 



dt / p''dx< 



Observe that ad/ 2 < 1. While inequalities of this kind are well known to the experts [17] . 
the references that include the case of exponents less than one are not common. For the sake 
of completeness, we provide a sketch of a simple proof of inequality (14. lip in Appendix II. 
Therefore, from (14.101) we can conclude that 

Now, suppose that C{q)x J^2 p{x, t) dx drops below 2 at some time to. Then, for all later times, 
we get 

dt [ p'^dx<-C{q) [ \Vp''^^\^dx (4.12) 

(we use C{q) for a positive constant depending only on q that may change from line to line). 
Let us recall another Gagliardo-Nirenberg inequality 



/ \Vp''/^\^dx(4-C{q)x [ pdx)-qe[ p^'^-Ux. 



II^^L/''-^' <C(g,rf)||Vt;|U2||t;||^<'-'. (4.13) 
Applying it in (I4.12p with v = p^^^ in = 2 leads to 

dt [ p'dx < -C{q) ( [ p'^dx^ ' ' ( [ pdx] " ' ■ 

Solving this differential inequality, and using the fact that J^^ P dx is monotone decreasing, 
leads to 

/ p{x, ty dx < min / p{x,toydx,C{q){t-to)~'^+^i / p{x,to)dx) ]. 
Then the argument identical to that in the proof of Theorem 11.11 implies that 

( 9-2 \ 
^l|p(-,to)|Ui,C(g)e-^ [j^r^^y^j (4-14) 

(observe that the proof of Lemma [2l2] goes through when C{q)x J^2 p{x,t) dx < 2). Since we 

have a uniform in time upper bound ||p(x, t)||Loo < max((x/e) ||po||l°°) (see Lemma [5l6] 
below), (14.141) implies the first statement of the theorem. Moreover, we can always take initial 
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data such that to = 0, and the L°° norm of po is sufficiently small (< (x/e)?^), making the 
bound on the right hand side of fl4.14p equal to c{q)x^^- This proves the second statement of 
the theorem. □ 

5. Appendix I: Global existence of smooth solutions 

Here we prove Theorem 14.11 We point out that equations involving chemotaxis and logistics 
or bistable type reactions have been considered by many authors (see e.g. [201 ESI ESI E21 EI] 
where further references can be found). In particular, in [211 EI], global regularity of solutions 
to a system similar to (II. ip was obtained. The main difference between these works and what 
we need here is that [291 !3T] work on a bounded domain. Since we consider a different setting 
and need different control of the solution (including moments) we present a brief sketch of the 
proof of global regularity result in this appendix. 

We begin with the construction of a local solution in an appropriate space. We will consider 
arbitrary dimension d. Recall that 

ii/iim„= / (ip(x)i + ivp(x)i)(i + |xnc?x, 

and the Banach space i^s,n is defined by the norm „ = ||/||Af„ + II/IIh"- First, we need a 

simple lemma on the heat semigroup action in this space. 

Lemma 5.1. Assume that po G -ft's.n, with s > 0,n > 0. Then we have 

l|e*^Po||M„ < C(l + r/^)||po||M„, ||Ve*^po||M„ < C{t~''' + t'~--'^l')\\p,\\M^- (5.1) 

We'^PoWus < ||Ve*^po||i/. < Ct-^'^Wp^H^. (5.2) 

As a consequence, 

l|e*^Pok.„ <C7(l + r/2)||po||^,,„, ||Ve*^po|k.,„ <C(t-V^ + t(--^)/^)||po||x.,„. (5.3) 
The proof of Lemma [5. II is elementary and we omit it. 

Next, we set up the contraction mapping argument for local existence. We will use the 
Banach space Xj„ = C{Ks^n, [0, T]) with a sufficiently small T > 0. Let us rewrite the equation 
(14. ip in an integral form using the Duhamel principle. 

p{x, t) = e*^po(x) + / e(*-^)^ (-V • (up) - ep^ + xV ■ (pVA-^p)) ds. (5.4) 
Jo 

Let us denote 

B,{p) = f e(*-^)^ (-V ■ {up) - ep^ + xV ■ (pVA-^p)) ds. 
Jo 

We need the following auxiliary estimates. 

Lemma 5.2. Assume that q, s,n are positive integers and s > | + 1. Let f,g E Kg^n- Then 

\\P - q'Whs < CiWfr^} + ||^?||?,-^)||/ - g\\H^ (5.5) 
||/VA~V - gV^^'gWus < cmUs + \\g\M\\f - gWu^ (5.6) 
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\\r - 91m. < chiWhs' + wgrns'm - 9\\m. (5.7) 

ll/VA^V - gVA-'g\\M„ < + ||^?||^,. + ||^7||mJ(||/ - 9\\m. + \\f - gM- (5.8) 

All constants in the inequalities may depend only on q, d, s and n. 

Proof. All these estimates are fairly straightforward. The estimate (15. 5p follows from writing 
ji — gi = (^j — g)[f''^^ + ■ ■ ■ + g'^^^) and the fact that if** is an algebra when s > d/2 (see, e.g. 
[35]). The estimate (15. 6p follows from a similar argument. The third inequality (15.71) is proved 
by the same expansion and use of Sobolev imbedding implying + ||V/||l°° < C*||/||_h''' 

and similar bounds for g. Finally, to prove the last inequality (15. 8p . write 

/VA- V - gVA-'g = (/ - g)VA-'f + gi\/A-'f - VA-'g). 

Integral of the right hand side expression against (1 + |x|") does not exceed 

||/-^7||AfJ|VA-VlU- + lkllA/J|VA-i(/-^7)|U- <C'(||/||h= + ||(7||mJ(||/-^7||m^^ 

For the case of the gradient, observe that 

V ■ (/VA-V - g\/A-'g) = (V/ ■ VA-^/ - \/g ■ \/A-'g) + (f - g'). 

The first two terms are then controlled similarly to the previous estimate, while the last two 
terms are easy to handle. □ 

Now we can prove a key Lemma setting up contraction mapping. 

Lemma 5.3. Suppose that u G C°°(]R'^ x [0,oo)) and V ■ m = 0. Let s,q and n he positive 
integers, s > | + 1. Let f,g G Xl,n- Then 

WBrif) - BTig)\\xT^ < "11/ - 9\\xt., (5.9) 

where for T < 1, we have 

a < C{d,q,n,e,x) max {\H;t)\\c^ + ll/(-,^)llC + U-M'C + ll/(-> + lk(-,^)k.,„) T^^'- 

(5.10) 

Proof. Consider 

B,{f) - Bt{g) = f e^(*-^) (V(«(/ - g)) - eif - g') + xV(/VA- V - g^A'^g)) dr. 
Jo 

Using Lemmas 15.11 and 15.21 we find 

\mf) - B,{g)\\K.^^ <C f [{{t - r)-V2 + (t _ (||^||^, + ||;||^^^^^ + ||^||^^^ J 

+ (1 + - rT") (ll/fxl + IkllC)] 11/ - ^11^^.^ 

< c [{e" + 1("+^)/2) maxo<.<, (ll^llc^ + + Ikk.J 

+ (t + t("+2)/2) maxo<.<, (ll/fx:! + ^^^o<r<t\\f - g\\K..„. (5.11) 

Therefore, we obtain (15. 9p . For T < 1 we can ignore the higher powers of T and the estimate 
for a follows from flSHD . □ 
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In a standard way, Lemma [5^ implies existence of local solution via the contraction mapping 
principle. 

Theorem 5.4. Assume q,s,n are positive integers and s > | + 1, -u G C°°(M'^ x [0,oo)), 
V ■ M = 0. Suppose po ^ Kg^n- Then there exists T = T{q,d,u, s,e,x, WpoWks.u) sttc/i that there 
exists a unique solution p{x,t) G Xj„ of the equation fl5.4p satisfying p{x,0) = po{x). 

Remark. Higher regularity of the solution in space and time (in particular implying p{x, t) G 
C{H'^,{0,T]) for every m > 0) follows from Theorem 15.41 and standard parabolic regularity 
estimates applied iteratively. 



Corollary 5.5. If under conditions of Theorem \5.4\ we prove global a-priori estimate on t)]!//^ 
and \\p{-,t)\\Mnj then the local solution can be extended globally to X'^^ with arbitrary T. 

Indeed, if there is a control on growth H'^ and M„ norms of the solution, we can just extend 
it by iterative application of local result as far as we want. To prove the bounds on M„ and 
j^s ]2orms of the solution, we first establish control of the L°° norm. 



Lemma 5.6. Assume that p{x,t) is the local solution guaranteed by Theorem \5.4\ Then 

||p(-,t)IU- < iVo = max ((x/e)^, ||po||l-) (5.12) 

for allO<t<T. 

Proof. Assume this is false, and there exists A^i > A^o and < ti < T such that we have 
||p(a^, ^i)||l°° = Ni for the first time (that is, for all x and <t <ti, \p{x,ti)\ < Ni). We claim 
that in this case there exists Xq such that p(xo,ti) = Ni. Indeed, the only alternative is that 
there exists a sequence Xk such that p{xk,ti) — )■ A^i as /c — ?• oo. If x^ has finite accumulation 
points, set one of them as xq. By continuity p(xo,ti) will be equal to Ni. Thus it remains to 
consider the case where — )■ oo and passing to a subsequence if necessary we can assume that 
unit balls around Xk, Bi{xk), are disjoint. By a version of Poincare inequality (see e.g. [32]), 
we have \\p -p\\l^iB^(x,)) < C\\p\\l.^B^{x,)y Since J2k WpWl^iB^ix,)) < C'(ti) < oo, we get that 



1 



But this is a contradiction with f^^^ |p(x)|(l + |x|") rfx < C(ti). 

Therefore, there exists Xq such that p(xo,ti) = Ni (we consider the case of a maximum; the 
case of minimum equal to —Ni is considered similarly). Then 

dtpixo, t) 1^^^^ = {u ■ V)p{xo, ti) + Ap{xo, ti) + xVp(a;o, ti) • VA-'^p{xo, ti) 
+Xpixo, tif - ep{xo, tiY < p(xo, tif{x - ep(a;o, ^i)""^)- 
By assumption on Ni, we see that dtp{xo,ti) < 0, contradiction with our choice of ti. □ 

Let us first prove an upper bound on the growth of the M„ norm of the solution. 
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Lemma 5.7. Assume that p{x,t) is the local solution guaranteed by Theorem 5.4 ■ Then on 
the interval of existence, we have the following bound for the growth of the Mn norm of the 
solution. 

\\pi;t)\U < C(l + r/')||po||M„exp (c f [e(l + r-'')\\p{.,r)\\l-^+ (5.13) 



dr) 



Proof. Consider (15.41) . By LemmaEHl we have ||e*^Po||M„ < C'(l + ||polUf„- Also, note that 

\\up\\m„ < \\u\\cA\p\\m„, Wp^Im^ < ||p|li^^||p|U/„, and 



||pVA-VIIm„ < c||va-VIU-IIpIIm„ < c{\\p\\l^ + \\p\\li] 

Therefore, applying these estimates and Lemma [5. II to (15. 4p . we obtain 

l|p(-,t)llM„<C(l + r/2)||po||M„+C r([e(l+r"/2)||p(-,r)||ri+ 

^0 

(,-1/2 + r(--')/^){\H;r)\\c. + ||p(-,r)|U^ + ||p(-,r)|UO] llp(-, OHmJ dr. 



(5.14) 



The inequality (15.141) implies the bound (15.131) . This upper bound is not optimal, but it is 
sufficient for our purpose. □ 

Now we are ready to prove uniform in time bounds on the H'^ norm of the solution. 



Lemma 5.8. Let p{x,t) be the local solution whose existence is guaranteed by Theorem 5.4 
Suppose that ||p(-, t)||L°° does not exceed Nq for all < t < T. Then 



\\p{-,t)\\Hs < max(||po||H=,C(u,d,g, s,x,e, A^o)) 



(5.15) 



Proof. Consider for simplicity the case where s is even (the odd case is very similar). Apply 
A^/^ to (14.11) . multiply by A^/^p(x, t) and integrate. We obtain 

^dMh = J [A^/'iu ■ V)p](A^/V) dx-ej (A^/V'')(A^/V) dx - \\p\\U^ 

+X [ [V ■ A'^''{pVA-^p)]{A'/^p) dx. 
Using V ■ M = 0, we obtain 

[A'^\{u-V)p)]{A'^^p)dx <C\\u\\c 



(5.16) 



Next, the second integral on the right hand side of (I5.16P can be written as a sum of a finite 
number of terms of the form J^^ -D^'pHLi D^^pdx, Si + ■ ■ ■ + Sg = s, Si > 0. Here denotes 
any partial derivative operator of the Ith order. By Holder's inequality, we have 



D'pYlD'ydx 



1=1 



i=l 
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Yli=iPi ^ = 1/2- Take pi = 2s/ Si, and recall the Gagliardo-Nirenberg inequality ([TOl |22l [T7] ^ 



\D'^p\\l2s,s^ < CIIpII^oo^ \\D'p\ 



L2- 



(5.17) 



Then we get 



<c||p|II 



«-l||^||2 

PWh^ 



Finally, we claim that the third integral on the right hand side of fl5.16p can be written as a 
sum of a finite number of terms of the form Jj^^D^pD^pD^^'^^'^A^^pdx, where k = 0, . . . ,s. 
The only term one gets from the direct differentiation that does not appear to be of this form 
is J^a A*/^pVA'*/^pVA~"^p(ix. However, integrating by parts, we find that this term is equal to 
-lJ^,\A'/^p\^pdx. Now 



D'pD^pD'+'^-^A-^pdx 



<c\\d'p\\lA\d'p\M\d^-^p\\l^. 



Pi^ + P2^ = 1/2, p2 < CO- Here we used boundedness of Riesz transforms on L^^, p2 < oo. Set 
Pi = y , P2 = ^3^- By Gagliardo-Nirenberg inequality f l5.17p with Si = k, s — k, we get 



[ D'pD'^pD'+'^-^A-^pdx 



< C\\p\ 



2 

Ha- 



Putting all the estimates into fl5.16p . we find that 



■^dt\\p\\]js < <^||p||l-||pIIh= 



I l|2 ^ r^w II II ||2 

\P\\h=+^ — IIPII-L°° IIPII//^' 



2+ 

\p\\hs 



-d/2 



-d/2 



(5.18) 



We used another Gagliardo-Nirenberg inequality in the last step. The differential inequality 
( I5.18P implies the result of the lemma. □ 

Given the M„ and norm bounds we proved, the solution can now be continued globally, 
completing the proof of Theorem 14.11 



6. Appendix II: The Gagliardo-Nirenberg inequality with p <1 

Twice in the paper, we needed to apply Gagliardo-Nirenberg inequalities with one of the 
summation exponents less than one (see (14. lip . (I4.13P ). Such inequalities are certainly known 
and can be found in mathematical literature (see e.g. encyclopedic [IZ]). However, it was 
not easy for us to find a reference with a transparent self-contained proof, and for the sake 
of completeness we provide a sketch of an elegant and simple proof here. The idea of this 
argument has been communicated to us by Fedor Nazarov. We will prove a slightly more 
general inequality containing both (14. lip and (I4.13p . 



Theorem 6.1. Let f G Cn 



l'^),d> 2. Then 
<C(g,d)||V^;||22||l;||i7^ 



1 _ 1 + 1- 

d 2 ^ r 



(6.1) 



The inequality holds for all q,r > such that q > r and ^ — + - > 0. 
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Proof. Let denote regions in M.'^ such that \Ak\ = 2^^^, G Z, the boundary of A^. coincides 
with a level set of \v{x)\ = ffc+i, and \v{x)\ > Vk+i inside A^. Then 

\\v\\h < "^IM^l 

fcez 

Fix a small S > 0. Let us call k "important" if v^+i < (1 — S)vk. Denote the set of all important 
k hj I. Observe that 



keZ k€l 

Indeed, a sequence of not important consequent k contributes at most ^;->o ^^'^'(l— 
compared to the contribution lAfc+ilt;^^^ of the single next term. 
For the U norm, we have the estimate 

{l-a)q/r 

^^'-^'>c{j2\A,K 

For the gradient term, by the co-area formula (see e.g. [5J) we have 

j \Vv\dx= V.'^-^{x: \v{x)\ = s)ds, 

where 'H'^~^ is the d — 1-dimensional Hausdorff measure. By the isoperimetric inequality, 

'W'-^x : \v{x)\ = s)> C\Ak\^--^ 
if s > Vk+i (see e.g. [5]). Therefore, 

/ l^^vldxyClAkl^-^Vk-Vk+i). 

Jvl-,1 <v(x)<vu 



\V\\Lr 



By Cauchy-Schwartz, 



f l^^vl^ dx > I f \Vv\dx] >C\Ak\^-^Vk-Vk+iy. 



Therefore, 



fcez kei 
Thus, it remains to prove that 

(\ {l-a)q/r / v aq/2 

k£Z / \kel / 



(6.2) 



Observe that, if > 3, then we have 



aq/2 ' - ^2(^ 



.kai ) \k&i 



2d 



2d 
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(since ^f^bl > {EkbkY for 6^ > 0, < s < 1). Write 



14, \ {l-a)q/r {l-a)q 
\^k\ Uk 



aq(d-2) 
2d 



(6.3) 



Apply Holder inequality on the left hand side of (16. 2 p , rasing the first term in (16. 3 p to the power 
^^^3^, and the second term to the power ^g^f_2) ■ Notice that the inverses of these powers sum 

to one due to the definition of a in (16. ip . The resulting inequality coincides with (16. 2p . Finally, 
when d = 2, we have a = 1 — q/r, and (16.21) follows from a more elementary consideration. □ 
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